function [CovBlocksNoisy] = addNoiseSamples(Samples, SampleLengths, CovBlocks, Combinations)

N=length(CovBlocks);

% We add for each block i symmetric Gaussian noise with variance i. 
CovBlocksNoisy = CovBlocks;
for k=1:N
    CovBlocksNoisy{k}=cov((CovBlocks{k}^0.5*Samples(Combinations{1,k},1:SampleLengths(1,k)))');
end